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Abstract 

Metallic foams are being investigated for possible use in the thermal protection systems of reusable 
launch vehicles. As a result, the performance of these materials needs to be characterized over a wide 
range of temperatures and pressures. In this paper a radiation/conduction model is presented for heat 
transfer in metallic foams. Candidates for the optimal transient experiment to determine the intrinsic 
properties of the model are found by two methods. First, an optimality criterion is used to find an 
experiment to find all of the parameters using one heating event. Second, a pair of heating events is 
used to determine the parameters in which one heating event is optimal for finding the parameters 
related to conduction, while the other heating event is optimal for finding the parameters associated 
with radiation. Simulated data containing random noise was analyzed to determine the parameters 
using both methods. In all cases the parameter estimates could be improved by analyzing a larger 
data record than suggested by the optimality criterion. 
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Nomenclature 


a — Coupling coefficient, m K /W 3 
b = Heater thickness, m 
c — Weighted specific heat, J/(kg K) 
c g = Specific heat of gas, J/(kg K) 
c/t = Specific heat of heater, J/(kg K) 
c m = Specific heat of metal, J/(kg K) 

Dj = Fiber diameter, m 
d g = Gas collision diameter, m 

e = eO 4- ep * T, Specific extinction coefficient, m 2 /kg 
F = Efficiency factor for conduction 
K b = Boltzmann constant 
Kn — Knudsen number 
k — Thermal conductivity, W/(m K) 
k g = Gas thermal conductivity, W/(m K) 
k* — Gas thermal conductivity at 1 atm, W/(m K) 
km = Thermal conductivity of metal, W/(m K) 

6 — Characteristic length, m 

P — Pressure, Pa 

Pr = Prandtl number 
T = Temperature, K 
T s = Sampled temperature, K 
t — Time, s 
x — Spatial coordinate, m 
a — Thermal accommodation coefficient 
(3 — Specific extinction coefficient, W/(m 3 K 4 ) 

e — Porosity 

6\ — Emittance of the septum plate 

£2 — Emittance of cold plate 

7 — Specific heat ratio 

A m = Molecular mean free path, m 
q c = Conduction heat flux, W /m 2 
q r — Radiation heat flux, W/m 2 
qr = Total heat flux, W /m 2 
p g = Density of the gas, kg/m 3 
Ph ~ Density of the heater, kg/m 3 
p m = Density of metal, kg/m 3 
Pf = Density of the foam, kg/m 3 
a — Stefan-Boltzmann constant 
uj = Albedo of scattering 
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1 Introduction 


High porosity metallic foams are being studied as possible components of aero-space thermal protection 
systems. The performance of these systems needs to be characterized over a wide range of temperatures 
and pressures during ascent and re-entry. 1 In the future when such materials are specified as part of a 
vehicle program, part of the procurement process will involve certification that the material meets the 
specifications. 

To date there has been little research on experiment design for thermal characterization of high-porosity 
materials. The present research is intended to close this gap in the procurement cycle by investigating 
transient methods for accurate measurement of thermal properties in metal-foam materials. In this paper 
only numerical simulations are presented. 

A review of the pertinent literature is given next in the areas of optimal experiment design and heat 
transfer in high-porosity materials. Transient experiments combined with parameter estimation have been 
used for obtaining thermal properties for many years. 2 In these methods the desired parameters are found 
by non-linear regression between the experimental data (temperatures in this case) and a computational 
model of the experiment. Parameter estimation concepts have recently been applied to optimal experiment 
design for thermal characterization of uniform materials 3 and for materials with temperature- varying 
properties.' 1 One of us (Cole) has previously studied optimal experiment design for measurement of thermal 
conductivity in low-conductivity materials 5 and in layered materials. 6 The work was limited to a small 
rise in temperature so that radiation heat transfer was negligible. An optimality criterion was used to find 
the best experimental conditions for estimation of thermal properties. 

There have been many studies of glass-fiber insulation in which combined radiation and conduction 
heat transfer are present. Yuen et. al 9 used a detailed radiation model and measured optical properties to 
simulate the temperature response of a fibrous insulation material from first principles. The model involves 
a detailed determination of wavelength-dependent radiation exchange between multiple zones through the 
thickness of the material. Their results compared favorably to transient experimental tests carried out in 
1974 on LI9000 sintered silica insulation (rigid space shuttle tile). 

Insulation models based on optical properties, although science-based, do not provide direct evidence 
of the insulation’s ability to withstand reentry heating. Simplified insulation models, combined with direct 
tests, can provide compelling evidence that a thermal-protection system is ready for use on a human-piloted 
vehicle. Daryabeigi 7 used steady experiments to characterize low- density fibrous insulation over a wide 
temperature and pressure range and a transient experiment was used to simulate re-entry aerodynamic 
heating conditions. A combined radiation /conduction model was used, with the radiation described by 
a two-flux method assuming anisotropic scattering and a gray medium. Two models were discussed for 
combining solid and gas conduction into an effective conductivity, a parallel model and a model based 
on fiber orientation. Thermal properties estimated with both models yielded similar results, so it was 
not possible to choose one conduction model over the other. Low Rayleigh numbers and experimental 
verification ruled out free-convection heat transfer in the pores. 

There arc few studies of open-cell metallic foams. Zhu et. al 10 carried out simulations of a titanium 
foam material for the purpose of finding the minimum weight of the thermal protection by varying the pore 
size of the foam across its thickness. Their steady-state model included effective conductivity and radiation 
parameters. Compared to uniform materials, a density-graded material provides the same protection with 
less weight, with greater percentage savings for thinner insulation. 

Sullins and Daryabeigi 11 studied a Nickel-foam material over a wide range of temperatures and pressures 
with steady-state experiments combined with parameter estimation methods. The model for the heat 
transfer included a two-flux model for the radiation and combined gas/solid conduction. The effective 
conductivity was described by parallel gas-solid conduction, temperature-dependent gas conduction, and 
an efficiency factor for solid conductivity (as a small fraction of the bulk- metal conductivity). A gas-solid 
coupling term was also added to take into account an observed increase in conductivity at higher pressures 
and temperatures. The measurement of the model parameters was carried out with a non-linear regression 
technique, and the parameters reported. 

In this paper the methods of optimal experiment design are applied to a high-porosity nickel foam 
material, with thermal property values taken from Sullins and Daryabeigi. 11 Both conduction and radiation 
heat transfer are included. Based on a one-dimensional model of transient heat transfer, a large number 
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of simulated experiments have been studied to determine which experimental conditions provide the best 
estimates of the thermal properties. No laboratory experiments are reported. 

The paper is divided into several sections, as follows: analytical model; numerical solution; simulated 
experiments; optimality criterion; optimal experiment design; estimation of parameters; and conclusions. 


2 Analytical Model 


In this section the heat transfer model of Sullins and Daryabeigi 11 ’ 12 is used to model metallic foams. 
It is a one dimensional, two- flux model taking into account both conduction and radiation heat transfer. 
Contained within the model are five parameters considered to be intrinsic properties of the material that 
must be determined empirically. Three of these intrinsic properties are related to the radiative heat transfer 
and two are related to heat transfer by conduction. 

The transient heat equation with radiation is given by 


P dt dx dx ’ 


( 1 ) 


where T is the temperature, q" is the conduction heat flux, q" is the radiant heat flux, p is the density, 
and c is the specific heat. The conduction heat flux is given by 



(2) 


where k is the thermal conductivity. 

The gradient of the radiant heat flux is given by 

^=0{l- u )(G-4aT*), (3) 

where G is the incident radiation, w is the albedo of scattering, 0 is the extinction coefficient, and a is the 
Stefan-Boltzmann constant. The albedo of scattering is to be determined experimentally. The incident 
radiation is related to the radiant heat flux by 


,,__±dG 
Qr 3(3 dx ' 


( 4 ) 


Thus the heat equation (1) becomes 


pc 


dT 

dt 


dx \ dx ) 30 dx 2 


The incident radiation, G , satisfies the second order differential equation 

1 d 2 G 


30 2 (l - to) dx 2 

on the interior. On the boundary, G must satisfy: 

2 


+ G — 4<tT a 


( 5 ) 


(6) 


30(A) 

2 

30(a) 



= 4aT£, x = 0, 

(7) 

an 


= 4oT*, x = L, 

( 8 ) 


where is the temperature of the heater, T c is the temperature of the cold plate (300.6 K), — 0.85 is 

the emittance of the septum plate, and €2 = 0.92 is the emittance of the cold plate. 
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Next the material properties will be discussed for a specific high porosity nickel foam in a N 2 atmo- 
sphere. The product pc is the weighted average of the values for the metal (nickel) and the gas (nitrogen). 
Explicitly, 


pc = Pm Cm ( 1 - e) + PgCgt (9) 

The density, p g , and specific heat, c 5 , of the gas are temperature dependent and given in the appendix. 
For nickel, p m is taken to be 8900 and Cm = 444.0. The porosity, e, is 0.968 for nickel foam. 

Heat conduction is assumed to occur in parallel in the gas and solid with some coupling in the form 

k = ek g - h (1 - e) Fk^ + a ( Fk g k m ) 2 (10) 


where F is the efficiency and a is the coupling coefficient. The thermal conductivity of the metal is a 
function of temperature and is interpolated from tabulated data. 13 For the gas 


k 


9 ~ 


k* 

Z’ 


( 11 ) 


where k* is the temperature dependent conductivity of the gas at atmospheric pressure (see Appendix) 
and Z is given by 


Z = $ + 2$ 


2(2 — a)^Kn 
0(7 + l)Pr 


( 12 ) 


In this expression, the thermal accommodation coefficient, a, is taken to be one, 7 = 1.4 is the specific 
heat ratio, Pr is the temperature dependent Prandtl number (see Appendix), and Kn is the Knudson 
number. The Knudsen number is given by 


Kn = ^ (13) 

0 

where A is the gas molecular mean free path and S is the pore size of the metallic foam. The quantities 4> 
and were defined 11 as 


•©1 

11 

CD m 

Kn < 10 
Kn > 10 

(14) 

* - {5 

Kn < .01 
Kn > .01 ' 

(15) 


This was to provide an approximation of Z in the different ranges by neglecting one term when it was much 
smaller than the other. However, determining the optimum experiment requires finding the derivative of 
the temperature at the sensor with respect to the empirical parameters. Having a discontinuous expression 
for the gas conductivity resulted in inaccurate results for these derivatives when the simulated experiments 
were done near the boundary between two of the ranges for Kn. Therefore in this work 4> = = 1. 

The pore size for the metallic foam was assumed to be in the same form as given by Daryabeigi 12 and 
is given by 
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7T Df 

4 


(16) 


In this expression Df = \AE — 5 m is the diameter of the strut. The mean free path is given by 

KbT 
sj 2nPd 2 g 


(17) 


in which Kb = 1.38£ — 23 is the Boltzmann constant, P is the pressure, and d g = 3.798£ — 10 m is the 
gas collision diameter. 

The extinction coefficient is given by 


P = e- p 


(18) 
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where p is the density of the foam and e is the specific extinction coefficient. The specific extinction 
coefficient is taken to be a linear function of the temperature. Explicitly, 

e = eo 4- e\ T (19) 

The quantities eo, and ei are considered to be intrinsic properties of the media to be determined experi- 
mentally. 

Thus, for this model, heat transfer in a metallic foam is determined by the intrinsic properties F, a, 
eo, ei and uj. The goal is to find the optimal transient experiment to determine these parameters. The 
variables considered for different experiments are the heating power, sensor location, heating time and 
pressure. 


3 Numerical Solution 


The heat equation (5) and the radiation equation (6) are solved numerically. For each time step of the heat 
equation, the radiation equation (6) is discretized and solved. The solution is then used in a discretized 
form of the heat equation (5) to determine the temperatures for the next time step. The temperatures are 
then replaced with the new temperatures and the process is repeated. 

Discretizing (6) through (8) gives the tri-diagonal system of equations 
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where G{ and 7* are the radiation and temperature, respectively, at node i. The elements of the coefficient 
matrix are given by 


Ai 

Ci 

Bi 


2 
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( 21 ) 

( 22 ) 

(23) 


In discretized form the heat equation (5) is 


p c j _ i 


Ax, 
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(24) 


3 Pj Ax j 


GJVO- 


where superscripts have been added to T and G to indicate the time step. Solving for TJ' +1 gives 
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The material properties of the metallic foam, pc, k , and 0, are all functions of temperature given by 
either polynomial approximations or interpolated from tabulated data. The temperature used for £* — piCi , 
ki , or 0i is the average temperature of the two endpoints of the interval Ax*. Explicitly for £ = £(T) 


Ci = C 


Ti + Ti 


i+i 


( 26 ) 


4 Simulated Experiments 


In this section the simulated experiment for determining the thermal properties of the metallic foam is 
described. The simulated experiment consists of heating the (one-dimensional) sample of nickel foam, at 
some pressure, for a particular time period with a heater of a specified power, while taking temperature 
measurements at one location in the body. Thus the variables of the experiment are the heater power, the 
sensor location, the pressure, the heating time, and the total duration of the experiment. 

Increasing the power to the heater increased the ability to determine the material properties. Since 
there would be physical limitations on the heater power supplied, rather than using the power of the heater 
as a variable between different experiments, an approximation of a heating apparatus used by Daryabeigi 14 
was made. This apparatus consisted of 2500 W cylindrical heating elements spaced 1.27 cm apart. This 
corresponds to about 310 MW/m 2 and an effective thickness of b — .00098 m. The heater and the sample 
were placed 5.08 cm apart. A cold plate, kept at room temperature, is considered to be along the other 
side of the sample. 

The temperature of the heater, node one in the discretized form of the heat equation, is determined 
from a lumped mass model. Under this model the heater must satisfy the differential equation 


PhChb 


dJ\ 

dt 


Qin Q.out i 


( 27 ) 


where p/ } , — 2600 kg/m 3 is the density of the heater and Ch = 800 J/(kg K) is the specific heat of the heater. 
The power supplied to the heater is Qi n — 310 MW/m 2 when the heater is powered and zero otherwise, while 
q oui is the heat loss to the sample and the surroundings. 

The losses from the heater include radiative transfer to the surroundings and to the sample. The 
radiative heat losses are given by 


<iout = e h a(Tjt - T 4 ) + €l a(T 4 - T 4 ), 


( 28 ) 


where th = .05 is the emittance of a water-cooled mirrored surface behind the heater, T 0 is the mirror 
temperature (300.6 K), and T s is the temperature of the surface of the sample. Gas conduction from the 
heater to the sample is also included, described by 

OT 

Qout = kg ■ ( 29 ) 

The total heat flow out of the heater is q out = q r out + q^uv 

A typical experiment starts at room temperature, then the heater is powered for a period of time while 
sampling the temperature at a particular location. Sampling continues for an additional period after the 
heater power is off. The heater model is incorporated into the transient finite difference code after the 
temperatures for each time step are found. 


5 Optimality Criterion 

The goal of this study is to examine a wide range of simulated experiments and determine the best 
experiment for determining the properties of the metallic foam. The suitability of individual experiments 
is determined from the sensitivity of the temperature with respect to the sought-after thermal properties 
of the sample, F, a, eo, ei and u j. Relabelling the parameters 6^, k = 1, . . . , 5, the normalized sensitivity 
coefficient for the kth. parameter, at the ith time step, is defined as 

X k (i) = b k ^±. ( 30 ) 
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Perturbing the parameters by a small amount, (1 + 8)b k , the sensitivity coefficients can be approximated 
by 


X k (i) a 


[Ti((l + 6)b k )-Ti(b k )] , 

- ,« 


1 S 


(31) 


The value of 5 — .0001 was found to give well-behaved values for X. 

In general an experiment is better if the sensitivity coefficients are larger. In addition, for experiments 
with more than one parameter the sensitivity coefficients must be linearly independent. For these reasons 
the optimality criterion used is 


D = 


|X T X| 

(nTl ax )p- 


(32) 


The sensitivity matrix, X, is defined by 


X T = 


' dr u dT , dT~ 
bl db 1 ,b2 db 1 ’"' ,b5 db 5 _ 


(33) 


where T is a vector of the temperatures at each sampling time. In (32), D is normalized by the maximum 
temperature, T max , squared and the number of sample temperatures, n , all taken to the power p (the 
number of parameters) in order to obtain a fair comparison among different experiments. 3 

The sensitivities are partial derivatives at a particular value of each parameter. The parameter values 
used for this study are given in Table 1, as reported by Sullins. 11 For all of the examples the thickness of 
nickel foam was 13.6 mm. 


Table 1: Parameter values for simulated experiments 


F a [(m K/W) 3 ] 

u e 0 [m 2 /kg] e 1 [m 2 /(kgK)] 

6.85E-3 3.89E2 

9.85 -2.63E-3 9.93E-1 


6 Optimal Experiment Design 

Our experience with thermal conductivity measurements, with no radiation present, suggested that it is 
best to heat the sample as rapidly as possible for a period of time while taking measurements at the 
surface closest to the heater. This is also the best approach for radiation/conduction materials. A variety 
of simulated experiments were explored by varying the pressure, the duration of the heating and the total 
duration of the experiment. Initially the simulated experiments were terminated at the time corresponding 
to the peak value of the optimality criterion, D. It was found, however, that extending the experiment 
duration increased the accuracy of the parameter estimation. For experiments at lower pressures radiative 
heat transfer is more dominant, while at higher pressures heat transfer due to conduction is more significant. 
It was thus expected that the optimum case would occur at moderate pressures since three of the parameters 
are related to radiation and two are related to conduction. 

The approximate conditions for the optimum experiment were found first by trial and error to find a 
rough value of the maximum of D, after which an optimization routine was run with the rough value as 
the initial guess. The maximum value of D found was 3.38E-16. This value resulted from an experiment 
conducted at 0.190 mm Hg in which the heater was run for 269 seconds and the temperature was sampled 
every 10 seconds for a total duration of 840 seconds. In Figure 1 the optimality criteria is graphed versus 
time for a few simulated experiments in which the optimum heating time was used and the pressure was 
varied around the optimal pressure. In Figure 2, D is graphed versus time with the pressure fixed at 0.19 
mm Hg for various heating times. As can be seen from the figure, the exact heating duration is not too 
critical. For heating times between 260 and 290 seconds the maximum value of D only varies by about four 
percent. In general if the sample is continuously heated the D value rises to a maximum and then begin to 
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decrease as the system reaches steady state. If instead the heater is shut off, at some time, the the value 
of D dramatically increases until the system again nears steady state (room temperature). The rise in the 
D values when the heater is shut off is so dramatic that in most cases it is impossible to see the rise in the 
D values during heating if the entire experiment is graphed on the same scale. In Figure 3 the same curves 
given in Figure 2 are repeated along with a curve corresponding to constant heating for 1200 seconds. 
The separation in the four curves prior to the heater being shut off is due to each experiment reaching a 
different maximum temperature which is included in the normalization. The best time for heating (269 
seconds) appears to be near the location where the D curve changes from concave up to concave down. In 
Table 2 the maximum value of D obtained for various pressures and heating times is summarized. 


Table 2: Maximum D values for experiments 


P [mmHg] 



Heating Duration [sec] 




250 

260 

269 

280 

290 

300 

0.09 

1.30E-16 

1.30E-16 

1.30E-16 

1.30E-16 

1.28E-16 

1.27E-16 

0.14 

2.69E-16 

2.71E-16 

2.72E-16 

2.71E-16 

2.69E-16 

2.67E-16 

0.19 

3.19E-16 

3.23E-16 

3.38E-16 

3.25E-16 

3.24E-16 

3.22E-16 

0.24 

4.32E-16 

4.11E-16 

3.93E-16 

3.78E-16 

3.65E-16 

3.53E-16 

0.29 

4.73E-16 

4.51E-16 

4.32E-16 

4.15E-16 

4.01E-16 

3.89E-16 


It should be noted that for the optimal experiment one would want to stop taking data as soon as D 
reached its maximum value. From Figure 1 or Figure 2 this is at 840 seconds. This time is also not too 
critical, since for any value between 830 seconds and 900 seconds D will still be within one percent of its 
maximum value. Thus for this model the optimal experiment for determining the material properties is 
for the sample to be heated for 270 seconds at a pressure of 0.19 mm Hg with sampling every 10 seconds 
until 840 seconds have elapsed. 

In Figure 4 the heat flux for the optimal experiment is given as a function of time. The total heat flux 
is approximately a step function while the radiative heat flux and heat flux due to conduction are of the 
same order, as expected. The temperature profile for the optimal experiment is given in Figure 5. One 
can see from the temperature profile the system has not reached steady state when the heater is shut off. 


7 Estimation of Parameters 

In this section the parameters are determined using three simulated experiments. Experiment 1 is the 
experiment found using the optimality criterion in the previous section. Experiment 2 is the same as 
experiment 1 except data is taken for 2000 seconds. Experiment 3 consists of two heating events concate- 
nated into one simulated experiment, where one of the heating events is dominated by conduction flux 
and the other is dominated by radiative flux. Data for all of the simulated experiments are obtained by 
calculating the exact temperatures using the values of the parameters given in Table 1 and then corrupt- 
ing them with random noise with a standard deviation of 0.08 K. The standard deviation was chosen to 
correspond with the maximum standard deviations reported by Sullins 11 for actual experiments in the 
same temperature ranges. A nonlinear Levnberg-Marquardt parameter estimation routine is used to fit 
the model parameters. 15 

For simulated experiment 1 the optimum experiment determined in the previous section was used. It 
consists of running the heater for 269 seconds and taking measurements until 840 seconds elapsed. The 
gas pressure is 0.19 mm Hg. The residuals between the exact and perturbed (noise added) temperatures, 
Tpert — Texact , and between the exact temperatures and those calculated using the estimated parameters, 
Test—Texact , are given in Figure 6. The residuals between the exact temperatures and those calculated with 
the estimated parameters is less than 0.5 K and is biased towards overestimating the exact temperatures. 
In addition the error in the calculated temperatures exceeds the random noise added to the exact values. 
The estimated values of the parameters along with their 90% confidence intervals are given in Table 3. 

Since at 840 seconds the temperature was still significantly above room temperature, for simulated 
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Table 3: Parameter values for simulated experiment 1 



eo 

ei 

D 

F 

a 

Exact 

9.850E+00 

-2.630E-03 

9.930E-01 

6.850E-03 

3.890E+02 

Estimated 

9.838E+00 

-2.619E-03 

9.930E-01 

6.582E-03 

4.696E-F02 

Percent Error 

-0.127 

-4.17E-01 

4.38E-03 

-3.91E00 

2.07E01 

90 % Confidence Interval 

6.551E-03 

4.922E-06 

3.152E-05 

7.956E-05 

2.821E+01 


experiment 2 the same heating time was used, but the temperature sampling was extended to 2000 seconds. 
By 2000 seconds the temperature at the sensor had cooled to 305 K. The results were much better both 
in terms of the residuals (See Figure 7) and the parameter estimation (see Table 4). The errors in the 
temperatures were less than .05 K while the parameters were found accurately within about one percent. 


Table 4: Parameter values for simulated experiment 2 



eo 

ei 

UJ 

F 

a 

Exact 

9.850E+00 

-2.630E-03 

9.930E-01 

6.850E-03 

3.890E-F02 

Estimated 

9.853E+00 

-2.633E-03 

9.930E-01 

6.874E-03 

3.819E+02 

Percent Error 

3.489E-02 

1.225E-01 

-1.233E-04 

3.521E-01 

-1.824E+00 

90 % Confidence Interval 

6.569E-03 

4.925E-06 

3.191E-05 

5.118E-05 

1.529E+01 


Table 4 shows that simulated experiment 2 estimated all of the parameters accurately, except for the 
coupling coefficient, a. Various other experiments with single heating events were tried, without success, in 
order to find the coupling coefficient more accurately. The coupling coefficient is only significant when the 
conduction through the gas is significant. This means a will be more important in experiments conducted 
at high pressure. In addition it should be easier to determine a from an experiment in which the radiative 
heat transfer is less important. Based on this information, simulated experiment 3 was constructed from 
a pair of heating events: one that was optimal for determining the conduction parameters; and one that 
was optimal for finding the radiation parameters. Noise was added to the data from this pair of heating 
events. The data sets were then concatenated as one data set for the parameter estimation. It should be 
noted that although optimal heating events were found for both the conduction and radiation parameters 
with the other parameters fixed, for the analysis of the concatenated data, all of the parameters were 
considered to be unknown. 

For the conduction-optimal heating event, increasing pressure was found to improve the estimation 
of a and F. However, an increase in pressure from 1 to 10 atmospheres increased the D values by less 
than one percent. Since it would be much easier to obtain a pressure of one atmosphere this was used 
for the conduction-optimum heating event. In this case, it was also found that using more power in the 
heater was better. Therefore the only variables to determine was how long to power the heater and how 
long to sample the temperature. From Figure 4 it can be seen that the flux due to conduction reaches 
a peak at about 20 seconds. At one atmosphere this peak moves to about 23.6 seconds, which was used 
for the conduction-optimal heating event. The temperature was sampled every second instead of every 
ten seconds. It was found that a data duration of 150 seconds provided good estimation of conduction 
parameters. Strict use of the D parameter would have resulted in a shorter experiment, but as with the 
single heating event the sample had not cooled very close to steady state by that time. 

For the radiation-optimal heating event, D is the largest for a vacuum (no gas conduction), Sullins 11 
reported using a minimum pressure of 0.0001 mm Hg so this value was chosen as a practically obtainable 
minimum pressure. For the experiments discussed so far, the initial heating had a significant conduction 
flux (see Figure 4). In Figure 8, the ratio of the conduction heat flux to the radiation heat flux is plotted 
for a few different heating times. It can be seen from the figure that if the sample is heated for 30 seconds 
this ratio is nearly zero for a large portion of the experiment. This corresponds to the experiment being 
dominated by radiation. For this reason the radiation-optimal experiment was to heat the sample for 30 
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seconds and then data was sampled for 100 seconds. Again a longer sampling time was used than the D 
value indicated was optimum. 

Simulated experiment 3 consisted of the conduction-optimal and radiation-optimal heating events com- 
bined together. The residuals for simulated experiment 3 are given in Figure 9. The estimated parameters 
and their 90% confidence intervals are given in Table 5. As can be seen in the figure and the table, the 
temperatures and the parameters are estimated accurately. 


Table 5: Parameter values for combined conduction and radiation experiments 



eo 

ei 

u> 

F 

a 

Exact 

9.850E+00 

-2.630E-03 

9.930E-01 

6.850E-03 

3.890E+02 

Estimated 

9.844E+00 

2.626E-03 

9.930E-01 

6.838E-03 

3.905E+02 

Percent Error 

6.271E-02 

1.666E-01 

-2.696E-05 

1.792E-01 

-3.960E-01 

90 % Confidence Interval 

1.944E-02 

1.668E-05 

2.818E-05 

9.107E-06 

8.771E-01 


The sensitivity coefficients (normalized by their maximum value) are plotted for simulated experiment 
1 in Figure 10. The sensitivity for a has nearly the same shape as the sensitivity for F. This indicates the 
sensitivities are almost linearly independent and is probably the reason it is hard to determine the value 
for a. A different choice of conduction parameters could probably be found that could be determined more 
accurately from simulated or actual experiments. 

Simulated experiment 3 provides more accurate results than either experiment 1 or experiment 2. For 
all the experiments, the strict use of the normalization of D by the number of samples taken tended to end 
the experiments before all useful data was taken. Unnormalized values of X T X values would indicate the 
sample should be heated to steady-state and then cooled back to room temperature, which would require a 
very long experiment. The best data duration appears to be somewhere between the duration determined 
by strict use of the maximum value of the normalized and unnormalized values of X T X. 

8 Conclusion 

For future use of high porosity materials in thermal protection systems for reusable launch vehicles, their 
performance needs to be characterized over a wide range of temperatures and pressures. The goal of 
this research is to develop general-purpose protocols for designing transient experiments to measure ther- 
mal properties of high porosity materials. As a first step, a metal foam material was studied with an 
experimentally validated thermal model. 

The D optimality criterion and consideration of the heat fluxes can lead to suitable experiments for 
determining the properties. Maximization of the optimality criterion D lead to an experiment containing 
one heating event that could be used to determine the radiation and conduction parameters. Strict use 
of the D parameter leads to experiments that were too short in the sense that more accurate parameter 
estimates could be obtained from simulated data sets of longer duration. By considering the model itself 
and in particular its conduction and radiation behavior, a set of two shorter heating events were found 
that provided more accurate parameter estimates. 

The radiation parameter u was easily found to a high degree of accuracy by all of the experiments. 
The coupling coefficient, a, was the hardest parameter to determine. The sensitivity of the of the coupling 
coefficient and the sensitivity of the solid conduction efficiency, F, are similar in shape which indicates 
they are close to linearly dependent. The coupling coefficient was an ad-hoc parameter added to Sulims’ 
model in order to obtain agreement with experimental data at high pressures. This suggests that there 
may be a different pair of conduction parameters that could adequately model the conduction process and 
provide more robust estimation than exhibited here. 
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Appendix 

The temperature dependent properties for nitrogen are from Daryabeigi. 12 

8.1 Temperature Dependent Nitrogen Properties 


p g = -1.734 E - 4 + 342.216/T (34) 

Cg = 1.083.545 - 0.328T + 6.949E - 4T 2 - 2.82 E - 7 T 3 (35) 

k g = 2.048E - 3 + 8.751F - 5T - 2.462E - 8T 2 (36) 

Pr = 0.854 - 7.085F - 4T 4- 9.008F - 7 T 2 - 3.207P - 10T 3 (37) 
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Figure 1: D for 269 seconds of heating at various pressures in mm Hg 
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Figure 2: D for P=0.19 mm Hg and various heating times 
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Figure 8: Ratio of conduction to radiation flux for different heating times in seconds. 



Figure 9: Temperature residuals for simulated experiment 3. 
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Figure 10: Sensitivities for conduction portion of simulated experiment 3 
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